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Abstract 

We present the full next-to-next-to-leading order (NNLO) coefficient func- 
tions for the polarized cross section dAa/dQ for the Drell-Yan process p+p — > 
+' X' . Here 'X' denotes any inclusive hadronic state and Q represents 
the invariant mass of the lepton pair. All QCD partonic subprocesses have 
been included provided the lepton pair is created by a virtual photon, which 
is a valid approximation for Q < 50 GeV. Unlike the differential distribu- 
tion w.r.t. transverse momentum the dominant subprocess for the integrated 
cross section is given by q + q — > 7* +' X' and its higher order corrections 
so that massive lepton pair production provides us with an excellent tool 
to measure the polarized anti-quark densities. Our calculations are carried 
out using the method of n-dimensional regularization by making a special 
choice for the 75 matrix. We give predictions for double longitudinal spin 
asymmetry measurements at the RHIC. 

PACS numbers: 12.38.Bx, 12.38. Qk, 13.85. Qk 
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1 Introduction 



With the advent of the RHIC at BNL we have a new facility to study the 
spin structure of the proton, (for a review on the potential of the RHIC see 
[1]), which supplements the existing polarized lepton-hadron machines. Po- 
larized proton-proton collisions with a very high luminosity and a maximum 
centre of mass energy of y/s = 500 GeV will provide us with many more de- 
tails about spin distributions than possible with the existing lepton-hadron 
machines, which give very little information about the polarized gluon and 
sea-quark parton densities. If we limit ourselves to the Drell-Yan process 
where the initial state hadrons are longitudinally polarized then the gluon 
initiated subprocess dominates the quark process in d 3 Aa / dQ 2 / dy / dpr at a 
transverse momentum px > Q/2 and fixed rapidity y [2], [3]. This result, 
however, depends on the parton densities involved. For a study of the cross 
section d 3 Acr / dQ 2 / dy / dp? and its dependence on the parton density set we 
refer to the complete NLO calculation in [4] (the non-singlet part was al- 
ready done in [5] and [6]). However for the totally integrated cross section 
dAa/dQ and for p + p — > l + l~ + l X' the sea quark initiated process domi- 
nates the reaction for the whole phase space. Therefore this process provides 
us with an excellent tool to measure the sea-quark densities. This process 
was calculated up to next-to-leading order (NLO) in [7] and subsequently 
confirmed in [8], [9] [10] and [11]. It is our goal to add all subprocesses 
which appear in the next-to-next-to-leading order (NNLO) contribution to 
dAa/dQ and to see whether the sea quark dominance persists. In order to 
calculate such a process one has first to specify the scheme for the 75 matrix 
if one wants to use n-dimensional regularization and massless quarks. In this 
paper we use the HVBM prescription of 't Hooft and Veltman [12] which is 
more elaborated upon by Breitenlohner and Maison in [13]. However this 
prescription like any others requires the introduction of evanescent counter 
terms before the reaction can be renormalized by the standard procedure. 
An example for the HVBM prescription is that the non-singlet axial vector 
operator gets renormalized in spite of the fact that it is conserved. This 
effect has to be undone by introducing an additional renormalization con- 
stant. This has also to be done for higher spin operators otherwise their 
anomalous dimensions do not equal the anomalous dimensions of the spin 
averaged operators. Also the Adler-Bardeen theorem [14] is violated which 
states that the Adler-Bell-Jackiw anomaly [15] does not receive higher order 
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QCD corrections. A scheme which is equivalent to HVBM is the one given 
by Akyeampong and Delbourgo [16]. Here the 75 matrix is replaced by the 
Levi-Civita tensor multiplied by four 7 matrices. If the Feynman and phase 
space integrals are calculated before the traces are taken the result is the 
same as in the complicated HVBM technique. More details are given in the 
next section. 

Our paper is organized as follows. In Section 2 we introduce our notations 
and discuss the technicalities which are involved when using Akyeampong- 
Delbourgo prescription for the 75-matrix. In Section 3 we present the NNLO 
corrections to the coefficient functions of the polarized DY cross section 
dAa/dQ. In Section 4 we study the NNLO corrections to polarized DY 
production in proton-proton collisions at the RHIC. Notice that this study is 
not complete since the three-loop contributions to the splitting functions are 
not available yet so that the corresponding polarized parton densities are not 
known. The long formulae for the NNLO coefficient functions can be found 
in Appendix A. 
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2 Kinematics of the polarized Drell-Yan pro- 
cess and 75 scheme 

The Drell-Yan process proceeds through the following reaction 

SO + H 2 (P 2 , S 2 ) -> 7*(?) +' , 

Z+(Zx) + Z"(Z 2 ) 

5=(P 1 + P 2 ) 2 , Q 2 = g 2 = (Z x + Z 2 ) 2 , (2.1) 

where Hi (i = 1, 2) represent the incoming polarized hadrons carrying the 
momenta Pi and spins Si. Further 1 X' denotes any inclusive hadronic state 
which is unpolarized. The lepton pair is represented by Z + Z~ with momenta 
Zi, Z 2 . In this paper we will only consider lepton pairs which have a sufficiently 
small invariant mass Q so that the photon dominates in the above reaction 
and Z-boson exchange effects can be neglected. The cross section is given by 



dAa^ 47ra 2 ATI , HlH „ ^ Q 
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A^ HlH2 (r,Q 2 ), r=^-. (2.2) 



dQ 2 3NQ 2 S v " S 



In the QCD improved parton model the hadronic DY structure function 
AW UlR2 is related to the coefficient functions A^ as follows 

AW^(r,Q 2 ) = E f-f ^A/f^^A/fO^ 2 ) 



T 



xA t3 ( — Q 2 lM 2 . (2.3) 
\x 1 x 2 J 

In the formula above Afi(x, jj 2 ) (i = q, q, g) are the polarized parton densities 
where fj, denotes the factorization/renormalization scale and x is the fraction 
of the hadron momentum carried by the parton. The DY partonic structure 
function AW ij is computed from the partonic subprocess 

i(Pi,si) +j(p 2 ,s 2 ) -> 7*(g) +ii(h) ■ ■ -i n (k n ) , (2.4) 

and it reads 

AW ij = K ij J d 4 q5(q 2 -Q 2 ) 
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i=l" V~"V ^-^i \ j=l 

x|AM t+ ^ r+n ...J 2 , (2.5) 

where denotes the colour and spin average factors and the polarized 
matrix elements are denoted by AM (when we refer to unpolarized struc- 
ture functions, matrix elements and parton densities we drop the A). Fi- 
nally note that the relation between the parton densities above and the 
parton momentum densities appearing in the parton density sets in the lit- 
erature or PDF libraries, which are denoted by A/ ( f DF (x, /i 2 ), is given by 
Af™ F (x,^)=x Af a (x,v 2 ). 

When lAMyl 2 in Eq. (2.5) is calculated up to order a 2 one encounters 
four partonic subprocesses which are characterised by the two partons in 
their initial state. In the case of quarks with a mass m ^ they are given 
by 

Q(Pi, s i) + Q(P2, s 2 ) -> 7* +' X' , 

\AM qq \ 2 = l - TV ( 75 ^ 2 fa - m)M l5 fa fa + m) M f ) , (2.6) 
9i(gi)(pi,si)+0(p2) ->7* +'*', 

\AM qg \ 2 = - e, uXa ^\ Tr (m" ^ fa ± m) M*) , (2.7) 
9i(9i)(Pi,«i) +q2(q2)(p2,s 2 ) ^7*+'X' 

| AM qiq2 \ 2 = I Tr ( 75 ^ 2 ± m) M 75 ^ (^ ± m) Aft) , (2.8) 
s(pi)+s(p 2 )^7* 



A2 7CT2 Al 7(Tl 

— e . P2 t2 e . Pl ^ TV (M^ 2 M u ^) (2 9) 



4 p 2 • h Pi- 1 
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where M denotes the matrix element which is given by the standard Feynman 
rules. Further the symbol TV can represent multiple traces when the matrix 
elements are calculated in higher order and for the reaction in Eq. (2.8) one 
must distinguish between qi = q 2 and qi ^ q 2 - The spin vectors Sj and the 
gauge vectors k (i = 1,2) satisfy the properties 

Si-pi = 0, Si-Si = -1, li-k = 0. (2.10) 

When the (anti-)quark is massless then one has to make the replacements 

7 5 ^(^i±m) -»• ±i 5 hifc, (2.11) 

where hi (i = 1,2) represent the helicities of the incoming (anti-)quarks and 
the + and — signs on the right-hand side hold for the quarks and anti-quarks 
respectively. The definitions above are chosen in such a way that the partonic 
polarized structure function satisfies the property 

AWij = Wij(+, +) - Wij(+, -) , (2.12) 

with +, — denoting the helicities of the incoming partons. 

The computation of the matrix elements in Eqs. (2.6)-(2.9) and their 
virtual corrections reveals divergences which occur when the momenta over 
which one integrates tend to zero (infrared), infinity (ultraviolet) or collinear 
to another momentum (collinear). The most popular way to regularize these 
singularities is to choose the method of n-dimensional regularization [12] in 
which the space is extended to n dimensions. The singularities are repre- 
sented by pole terms of the type (l/s) k with n = 4 + e. This method is very 
useful because it preserves the Ward identities in the case of gauge theories. 
However this is no longer the case when the 75 matrix and the Levi-Civita 
tensor appear like in Eqs. (2.6)-(2.9) or in weak interactions. There is no 
consistent way to generalize these two quantities in n dimensions contrary 
to the ordinary matrix 7^ or the metric tensor g^ v . In the literature one 
has proposed various methods to extend the 75 matrix and the Levi-Civita 
tensor to n dimensions but one always needs so-called evanescent counter 
terms [17] to restore the Ward identities. A very popular prescription is the 
HVBM-scheme which was proposed in [12] and generalized in [13]. In this 
approach the n-dimensional gamma-matrices and momenta have to be split 
into 4 and n — 4 dimensional parts. Therefore also the integrals over the final 
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state momenta have to split up in the same way. Many NLO calculations 
have been done in this scheme (see e.g. [5]- [11]). However this approach 
requires a special procedure to deal with the gamma-matrix algebra which 
is not implemented in the program FORM [18]. Since this program is used 
in our calculations we prefer another prescription for the 75-matrix which is 
given in [16]. It gives the same results as the HVBM-scheme but it is much 
simpler to use in algebraic manipulation programs. Moreover one does not 
have to split up the integrals over the final state momenta and one can sim- 
ply use the phase space integrals computed for unpolarized reactions. The 
procedure in [16] is given by 

1. Replace the 75-matrix by 

% % 
7 M 75 = g ^ paT l p l a 7 T or 75 = — e paTK Y Y Y 7* . 

2. Compute all matrix elements in n dimensions. 

3. Evaluate all Feynman integrals and phase space integrals in n-dimensions. 

4. Contract the Levi-Civita tensors in four dimensions after the Feynman 
integrals and phase space integrals are carried out. 



(2.13) 

Note that the contraction in four dimensions only applies to Lorentz indices 
which are present in the Levi-Civita tensors. The last step in (2.13) requires 
that one first has to apply tensorial reduction to all integrals. For simple 
expressions this takes the form 

r d n k 

J ^^*7(fc 1 Pi.ft) = A 00 (n)g^ + A 11 (n)p?p» + A 22 (n)p 2 l p 2 ' 

+A 12 (n)p1p 2 + A 2i (n) P2 Pi , 
JdPSWWk"f(k,p 1 ,p 2 ) = B 00 (n)g^ + B 11 (n)p?p» + B 22 (n)p 2 l p 2 ' 

+B 12 (n) j% p 2 + .821 (n) p% p\ , 

(2.14) 
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where the coefficients Aij, B, L j depend on n = 4 + e. However there exists 
an alternative method. If the matrix element only contains one 75 like in 
qg or even none as in gg one can apply a different rule. One replaces 75 
according to 1 in Eq. (2.13) and one contracts the Levi-Civita tensors in 
n dimensions before the Feynman integrals and phase space integrals are 
carried out. Further one factor of 1/2 in Eqs. (2.7) and (2.9) is replaced 
by l/(n — 2)(n — 3). The result is the same as given by the method in Eq. 
(2-13). 

The HVBM-scheme or the prescription above automatically reproduces 
the Adler-Bell-Jackiw anomaly (ABJ) [15] but it violates the Ward identity 
for the non-singlet axial vector current and the Adler-Bardeen theorem [14]. 
In order to obtain the correct renormalized quantities one therefore needs to 
invoke additional counter terms [19], which are called evanescent [17] since 
they do not occur for n — 4. This procedure was used to obtain the NLO 
anomalous dimensions for the spin operators which determine the evolution 
of the parton spin densities see [20], [21]. Evanescent counter terms are also 
needed in cases where collinear divergences show up like in partonic cross 
sections (see e.g. [22]). This is characteristic for the HVBM scheme as well 
as any other approach. The reason is that in n-dimensional regularization 
there exists a one-to-one correspondence between the ultraviolet divergences 
occurring in partonic operator matrix elements and the collinear divergences 
appearing in partonic cross sections provided both quantities are of twist 
two type [22]. In [23] we have calculated all evanescent counter terms for 
the HVBM-scheme up to second order. For later convenience we will express 
them in the unrenormalized coupling constant 



7 5,NS,- 



8(1 - x) + a s S £ 



'91 



e/2 



z {1) + ez {l) 

A qq ' qq 



+ aiSi 



'91 



47T 



(2.15) 



V5,NS- 
11 



'Q 2 



y (2),NS,- 
"11 



(2.16) 
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7 5,PS 



- n 2 9 2 \ — 



2\ £ 



y (2),PS 
"99 



(2.17) 



where /?o represents the lowest order coefficient in the beta-function. Note 
that we have to take lnQ 2 //^ 2 terms into account otherwise it is impossible 
to relate AW^ S to the unpolarized expression W^ s . This follows from the 
definition of 



■qq 



1 qq J w^(a s ,qy^) 



(2.18) 



Explicitly we can write 



5(1 -x) +a s S £ 



2\ e/2 



-8(1 -x)+e{ 



-8(1 - a;) ln(l - a;) + 4(1 - a:) In a: + 4 - 4x} 

+a 2 s S 2 e (^Pj C 2 F {- 16(1 -x) -(16 + 8x) In a: 
+ 16(1 — x) lnx ln(l — x)} 

+C A C F {i(^(l - x)) - - x) + 8(1 - x)C(2) 

+ ^ (1 - x) ln(l - x) + (-^ + 32 x) lnx - 4(1 - x) In 2 x} 



+n f C F T f { - - xj) + ^(1 - x) + 16 (1 - x) In 



x 



64 ^ 
ln(l-x)} 



(2.19) 



v5,NS,- 
99 



(j^J {C 2 F - l -C A C F ) 8(l + x)(4Li 2 (-x) 
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+4 \nx ln(l + x) + 2((2) - In 2 x - 3 In x) - 56(1 - x) 



(2.20) 

16(l-x)+8(3-a;)lnx + 4(2 + x)ln 2 x . 

(2.21) 

A comparison with previous calculations [22], [23] shows that the evanescent 
counter terms are universal except for the Q 2 / fi 2 terms and the e part of the 
order a s contribution. 



= alS 2 m n f C F T f 
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3 Corrections up to NNLO to polarized Drell- 
Yan production 

In this section we treat the various subprocesses which contribute to the 
Drell-Yan reaction up to NNLO. The Born reaction is given by the process 

q + (3.1) 

The first order corrections are given by the following subprocesses. The 
bremsstrahlung reaction to the Born process 

q + q^l* + g, (3.2) 

plus the one-loop correction to the Born process. Further we have a new 
process with a gluon in the initial state 

g + q(q) -»• 7* + q(q) • (3.3) 

In the next order we encounter the following processes. First we have the 
bremsstrahlung correction to reaction (3.2) 

q + q-*l* +g + g, (3.4) 

plus the two-loop correction to the Born process and the one-loop corrections 
to process (3.2). Next we have the bremsstrahlung correction to reaction (3.3) 

g + q(q) -»■ 7* + q(q) +g, (3.5) 

plus the one-loop correction to reaction (3.3). Subsequently we have three 
new reactions. The first involves a quark- ant i- quark reaction see Figs. 1, 2 

q + q^l* + q + q, (3.6) 

and the second contains two (anti-) quarks in the initial state see Figs. 2, 3 

q(q) + q(q) - 7 * + q(q) + q(q) , (3.7) 

In the above reaction the quarks can be equal or unequal. Finally we have 
a reaction with two gluons in the initial state 

g + g^l* + q + q- (3.8) 
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Figure 1: Annihilation graphs contributing to the subprocess q + q — » 7* + 
q + q. 



All matrix elements are calculated in n-dimensions using the algebraic manip- 
ulation program FORM [18]. The phase space integration and the Feynman 
integration are explained in [24] and [25]. The results can be expressed in 
the renormalization group coefficients APy, f3 and as follows 



-AW NS = 

qq 



5(1 - z) +a s S £ 



( 




AP4°)+<- ) +2,« + 



(91Y [1( 2 AP(°)®AP(°) 
l;U 2 / £ 2 \ qq w qq 



-2 ft AP(J) } + \ { AP«< NS <+ + 2 AP (o) W W _ 2 ^ W( 
+4 4J) ® APW - 2 ft zg> } + 2 AP £> ® - 2 ft «># 



99 



(2),NS 



+ 4 *£> ® AP W + 2 zg> ® «,« + 4J) ® z\ 



(1) 
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Figure 2: Gluon exchange graphs contributing to the subprocess q + q 
7* + q + q and q + q-^Y + Q + Q- 
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Figure 3: Gluon exchange graphs contributing to the subprocess q + q 
7* + q + q with identical quarks in the initial and/or final state. 



The qq reaction (see Figs. 1,2) yields 

— ApW)®ApW) + -t-ApW 

2 £ 2 99 w 99 ~ p \ 9 99 



AW, 



PS 
99 



PS 



+2 APW ® «,£>} + 2 APW ® tD« + ^' PS + 2 ,(f ps 



(3.11) 

The qq can be split into qiq2 (qi ^ 52) and a gg reaction. The latter does 
include Fig. IB. 



w 



(2)> p S _ „„(2) . (2) . (2) 



9192 



W CC + W DD + W CD , 



(3.12) 



w, 



(2),PS _ (2) , (2) , (2) , (2) , (2) 



(2) 



99 



W B B + W BC + W BD + WCC + + W CD 



(3.13) 
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Processes with unequal quarks in the final state (see Fig. 2) yields 



AW, 



9192 



2 £ 2 19 91 ~ P \ 9 «<? 



£ >-2 



+2 APW <} + 2 APW + + 2 ^ PS 



(3.14) 



to 



(2) 



(2) , (2) . (2) 
WCC + ™DD + <D 



(3.15) 

Processes with equal quarks in the final state (see Figs. 2, 3) yields 



AW qq = -a 2 S 2 ( % 



2\ e 



2 e 2 19 ^ 91 p I 9 99 99 



£ 1-2 



+2 AP(J) ® < } + 2 APW ® tDg) + < + 2 ^ PS 



_9 -(2),NS, 
ZZ 99 



(3.16) 



(2) (2) , (2) . (2) . (2) . (2) . (2) . (2) . (2) . (2) 

gq = ™CC + W DD+WcD + W FF + WeE + W FE + W C F + W CE + WdF 

(3.17) 



W 2) 



Finally we have the gluon-gluon scattering process 



AW, 



mi 



AF<«> 8 APi»> + - AP<„°> ® W <!> 



19 ^ 19 19 



+2AP(Sw (1) + w (2) 

1 93 ^ 99 1 99 



(3.18) 



where ® denotes the convolution symbol defined by 

f®g(z)=[ dzi [ dz 2 5(z - z 1 z 2 ) f{z{) g{z 2 ) 
Jo Jo 



(3.19) 
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The expressions follow from the renormalization group equations. They are 
constructed in such a way that they become finite after coupling constant 
renormalization and mass factorization are carried out. The (3q appears as the 
lowest order coefficient in the beta- function. The splitting function APij(z) 

and the coefficients wf-(z) with z = Q 2 /s also occur in the coefficient func- 
tions given below except for the NLO terms wff, which are proportional to 
e. They are given by 



w, 



(i) 
qq 



C F 8 V 2 {x) - 6 C(2) V (x) + (1 + x) ( - 4 ln 2 (l -x) + 3 C(2)) 

+ I±f! (_ 4 \ nx i n (i _ x ) + i n 2 x \ + 4 (x _ a-) + _ x ) ( 16 
1 — x v 



21 



~C(2)) 



(3.20) 



(i) 

?9 



i(^-l)(^(^)-3C(2 



x In 



x 



- 5 + Ax + x 



(3.21) 



where T>i denotes the distribution 



1-z 



(3.22) 



To render the partonic cross sections finite one has first to perform cou- 
pling constant renormalization. This is done by replacing the bare coupling 
constant by the renormalized one i.e. 



a s = a s (n 2 ) 



1 + a s (/i 2 ) S £ - p + 

e 



(3.23) 



To remove the evanescent counter terms we perform the following operation. 
We multiply each quark line with the inverse of Z^ q [22]. This means that 



16 



qq, qq and gi^2 are multiplied by (Z q ) 2 . The qg is multiplied by (Z t 



-qqj 



and the gg subprocess gets no factor. The remaining divergences, which are 
of collinear origin, are removed by mass factorization 

1 Q 2 



k,l=q,q,9 



r«l- 



A H (J) , (3.24) 



where ct is the number of (anti)-quarks in the initial state. Further Tki(z) de- 
note the kernels containing the splitting function which multiply the collinear 
divergences represented by the pole terms l/e k . For the different subprocesses 
the mass factorization relations in Eq.(3.24) become equal to 

r zL) 2 aw™ = rj? ® rg> ® a$> , (3.25) 



qq 



qq 



qq 



qq ' 



(z 5 qq ) 1 Aw qg = r qq ® v qg ® A qq + r qq ® v gg ® A qg , (3.26) 

(Z 5 qq ) ~ 2 AW qq = T qq ® T qq ® A 99 - + ® r gg ® A OT + T gq ® r 59 - ® A OT -(3.27) 



Kf 2 A ^< 



91 92 



+r 



992 



- 91 92 
r 



a + r 

^99 ~ 1 991 



9292 



A 992 



9J.71 ® 991 ' 



(3.28) 



99 



Tgg ® Agg 



Aqq + T 



.99 



+ r ? q ® T 9 q ® Aqq + Tqq 



Tqq ® Aqq , 



(3.29) 



AW, 



<-!<-! 



2T 



qg 



qq 



+T gg ® T gg 



Aqq + 2 Tqq 

) A n 



r^®Aq 9 +2rqq 



99 



*99 • 



■^99 
(3.30) 



Since we need the finite expressions up to order a 2 s it is sufficient to expand 
the kernels Y ki up to the following order in the renormalized coupling constant 

1 



■pNS . . -pNS 
99 99 



6(1 - z) + a s S £ 



-2ftAP£>) + lAP£> 



AP (o) 

qq 



NS,- 



— fAP(°)® AP(°) 



(3.31) 
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99 



Tgg CL S S £ 



2e q9 



+ a 2 s 



WW®APW-2A,APW) + 1(AP(J) 



(3.32) 



r PS = r P8 = r p ^ = r PS = i? 2 s* 2 

qq qq qq s 



— APW <g> AP(°) + — AP«- PS 



,(3.33) 



-nN_S _ pNS 
99 99 



-a 2 S 2 



J_ ad(1),NS, 

2 £ qq 



(3.34) 



gq ~ gq ~ u s 



a, S E 



-AP(°) 



99 



(3.35) 



r s9 = 5(1 -z)+a s S £ 



-AP(°) 



.9.9 



(3.36) 



After renormalization and mass factorization the coefficient functions have 
the following algebraic form 



N_S 
99 



-5(1 - z) - a s (n 2 ) 



APi°) In % + 



w 



(i) 



99 " ,,2 1 ^99 



{\ AP#> ® APW - \ (3 AP(J)} In 2 ^ + {AP«< NS 



-^o^+AP^O^jln 



Q 2 



2 ■ W 99- 



/i- 

(2),NS 



(3.37) 



A 99 



-AP(°) m^ W 1} 

^ 99 ^2 T 99 



(— AP(°)(8) f3AP(°) 
1 16 w V 99 
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(1) 

2 ~99 



+ Ap(f) - I ft AP«»} l„ 2 ^ + {\ap$ - ft <> + i „. 



i p - s 

qq 



(2) 
39 



(3.38) 



-a(/0 



{^>Ap(;)}ln 2 ^ + {iAP« 



PS 



Q 2 



+AP, ( g 0) ®< ) }ln^+ W , 



(2),PS 

qq 



(3.39) 



»9l92 



{^Ap(>Ap(;)}ln 2 ^ + {iAp(^ 



+AP, ( 9 0) ®< ) }ln 



Q 2 
/i 2 



w 



(2) 
91 92 



(3.40) 



A M = ~a 2 s (fi 2 ) 



(-AP(°) <g> APW) In 2 ^ + (- AP«< PS - AP «' NS '- 

\ g 99 ^ 99 J ^2 T \ 2 99 99 
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A 99 
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99 



(3.42) 



In Appendix A we give the explicit expressions for the coefficient functions so 



that one can determine the coefficients AP,^ and w£'. One of the features 



,(*) 
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is that the non-singlet part of the structure function satisfies the relation 



( z 5,NS,+\ 2 A ^NS = _^ 

V 11 ) m i 



NS 
<?<? 



(3.43) 



For the O(a^) correction this holds for Eqs. (A.1)-(A.13). In the presentation 
of the coefficient functions above we have put the renormalization scale fi r 
equal to the mass factorization scale /i. If one wants to distinguish between 
both scales one can make the simple substitution 



a s (/i 2 ) = a 8 ((j%) 



1 + 



47T 



■A) in 



t4 



(3.44) 
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4 dAa/dQ for the process p + p — > 7* +' X' 



In this section we will present total cross sections (see Eq. (2.2)) for polarized 
Drell-Yan production in proton-proton collisions at the RHIC and make a 
comparison with similar results in previous work. The cross section can be 
rewritten as 



In particular we study the dependence of the cross section on the input 
parameters such as the renormalization/factorization scale /i, the virtuality 
of the photon Q and the input parton densities. At this moment LO and 
NLO polarized parton densities are available but NNLO are not. Therefore 
an exact NNLO polarized cross section cannot be determined yet. Even 
an approximated polarized cross section cannot be given because a finite 
moment analysis of the anomalous dimensions is not yet available. Only 
the non-singlet anomalous dimension is known at least up to N = 14 [26]. 
Because of a lack of data the polarized parton density sets show a larger 
variation than in the unpolarized density sets. This in particular holds for 
the sea-quark and the gluon densities. For this reason polarized Drell-Yan 
production is very useful. The sea-quark contribution is dominant for the 
totally integrated cross section in proton-proton collisions while the gluon 
contribution dominates the differential distribution w.r.t. pr for px > Q/2. 
We shall use BB set 1 (BB1) for our plots [27] in the whole paper except 
for the longitudinal asymmetry where we shall choose also other density sets 
like BB set 2 (BB2) [27], the GRSV01 (standard scenario) and GRSV01 
(valence scenario) [28]. The gluon in BB1 is larger than in BB2. In its turn 
the gluon in the standard scenario is smaller than in BB2 but larger than 
in the valence scenario. For all these sets LO and NLO versions exist. It 
is clear that in the case of LO we will use the one-loop parametrization for 
the running coupling constant and for NLO the two-loop corrected running 
coupling constant. However for an estimate of the NNLO corrections we 




(4.1) 



where r = Q 2 / S and $y is the parton-parton flux defined by 




(4.2) 
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GRSV01 (LO, standard scenario) 


A^° = 


= 175 MeV 


S \ n J 


= 0.121 


GRSV01 (NLO, standard scenario) 


A NLO = 


= 257 MeV 


a^ LO (M z ) 


= 0.109 


GRSV01 (LO, valence scenario) 


A LO = 


= 175 MeV 


a^°(M z ) 


= 0.121 


GRSV01 (NLO, valence scenario) 


A NLO : 


= 257 MeV 


a* LO (M z ) 


= 0.109 


BB (LO, scenario 1) 


A LO = 


= 203 MeV 


a™{M z ) 


= 0.123 


BB (NLO, scenario 1) 


A NLO : 


= 235 MeV 


a^ LO (M z ) 


= 0.107 


BB (LO, scenario 2) 


A LO = 


= 195 MeV 


o%°(M z ) 


= 0.123 


BB (NLO, scenario 2) 


A NLO 


= 240 MeV 


a^ LO (M z ) 


= 0.107 


MRST02(LO, lo2002.dat) 


A LO = 


= 220 MeV 


o%°(M z ) 


= 0.125 


MRST01(NLO, alfll9.dat) 


A NLO 


= 323 MeV 


^ LO (M z ) 


= 0.113 


MRST02(NNLO, vnvalll55.dat) 


A NLO r 


= 235 MeV 


c£ NLO (M z 


) = 0.109 



Table 1: Polarized and unpolarized parton density sets with the values for 
the QCD scale A 4 and the running coupling a s (M z ). 

will only choose NLO polarized parton densities with the two-loop running 
coupling constant. The details are given in Table 1. Further we put rif = 4 in 
the coefficient functions and the running coupling constants and the densities 
are only presented for the u,d,s and g partons. Finally we also make a 
comparison with unpolarized Drell-Yan production. For this cross section 
we adopt for LO the MRST02 [29] densities, for NLO the MRST01 [30] 
densities and for NNLO the approximate MRST02 [29] densities, with one, 
two and three loop running coupling constants respectively. 

In Fig. 4 we have plotted the polarized Drell-Yan cross section in Eq. (4.1) 
up to NLO. It is clear that the qq contribution dominates the qg contribution. 
The former is negative above Q = 25 GeV, while the latter is already negative 
in a region above Q—13 GeV. This implies that the total NLO contribution 
is actually negative above Q=27 GeV where it is very small. For this reason 
we have plotted the absolute values of the contributions in Fig. 4. 

The relative sizes are not changed if we go to NNLO as demonstrated 
in Fig. 5. Again the qq contribution is negative above Q— 25 GeV and the 
qg contribution is negative between Q— 4 GeV and Q— 13 GeV. Here the 
qq and gg channels appear for the first time and they are both negative. 
However their contribution is negligible compared with the one given by qg 
and certainly by the qq result. Again we have plotted the absolute values 
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of the contributions. We conclude that via the qq subprocess the sea-quark 
contribution dominates the whole cross section in proton-proton collisions 
irrespective of the value taken by the gluon density. 

In Fig. 6 we have plotted the absolute values of the LO, NLO and NNLO 
polarized cross sections for 2 < Q < 30 GeV. Notice that we are in the range 
of small Q-values so the ^-factors are pretty large. We have defined the 
.fT-factors as 

a NLO Arr NNLO 

Aa LO ' A Aa LO ' ^ 

and plotted them in Fig. 7. At Q = 7 GeV the K-factors reach a minimum, 
j^nlo ^ and i^ NNLO ~ 1.3, and at larger Q they rise again. The rapid rise 
above Q— 18 GeV is due to the change in sign of the LO process near Q=25 
GeV where the ratio is infinite. However the magnitude of the polarized cross 
section is extremely small in this region. 

Next we turn to the variation of the cross section w.r.t. scale fi. In Figs. 
8a, 8b and 8c we show the LO, NLO and NNLO polarized cross sections at 
the scales fj, = Q/2, \x = Q and \x = 2Q. In each figure we see that if \x 
gets larger the cross sections increase at small Q but decrease it at larger Q. 
Therefore there are specific values in Q where the scale variations are small 
and as we go from LO to NLO to NNLO these points are at larger values of 
Q. The scale variation gets smaller as we go from LO to NLO as expected. 
In spite of the fact that we do not have the NNLO parton densities there is 
still an improvement in the scale dependence if we go from NLO to NNLO. 
The corresponding plots for the unpolarized cross section are shown in Figs. 
9a, 9b and 9c. These plots also show that when \x becomes large the cross 
sections increase at small Q and decrease it at larger Q. Further we see the 
expected overall decrease in scale variation as we go from LO to NLO and 
then to NNLO. 

To show the scale variation from a different point of view we plot the 
following quantity 

= ^4- (44) 

In Fig. 10a we plot the polarized quantity in Eq. (4.4) for 0.4 < /i//i < 2 at 
A*o — Q — 5 GeV. In this figure we see an improvement in the scale variation 
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while going from LO to NLO and then to NNLO. This feature also persists 
at higher Q-values i.e. Q — 10 GeV in Fig. 10b and Q = 15 GeV in Fig. 10c. 
For the unpolarized cross section the trend is the same. Here we have the 
NNLO parton densities and in Fig. 11a at Q = 5 GeV we see an improvement 
in scale variation while going to higher order. The same observation can be 
made in Fig. lib at Q — 10 GeV. However for Q — 15 GeV in Fig. 11c 
the NLO curve becomes slightly worse than the LO one. While on the other 
hand the NNLO curve is slightly better than the LO at least for larger scales 

Now we look at effect of the higher order corrections on the longitudinal 
asymmetry defined by 

All — ■ (4.5) 
a 

In Fig. 12 we have plotted All in percent for BB1 in the range 2 < Q < 30 
GeV in LO, NLO and NNLO. Since the polarized cross section changes sign 
we see that the asymmetry is negative for large Q (note the displaced zero). 
It is obvious that the longitudinal asymmetry is small reaching about one 
percent between Q = 10 GeV and Q = 20 GeV. 

To show the predictions of the various polarized parton density sets we 
compare the BB1 NLO contribution with the other three NLO parton den- 
sity sets i.e. BB2, GRSV01 (standard scenario or SS) and GRSV01 (valence 
scenario or VS) in Fig 13. Obviously the first three sets have small posi- 
tive longitudinal asymmetries while the GRSV01 valence scenario set has a 
large negative asymmetry. Therefore the set with the smallest gluon, here 
GRSV01 (valence scenario), has the largest asymmetry. On the other hand 
the set with the largest gluon, here BB1, has the smallest asymmetry. This 
feature is unchanged in NNLO as can be seen in Fig. 14. Hopefully the 
RHIC experiments will get sufficient luminosities and large enough proton 
polarizations to distinguish between the first three predictions and the last 
one. 

To summarize we have computed the NNLO corrections to the produc- 
tion of massive lepton-pairs in polarized proton-proton collisions. Hence as 
far as coefficient functions as concerned we have the NNLO corrections to 
both polarized and unpolarized processes. Unfortunately the corresponding 
anomalous dimensions of the partons are not available at NNLO. While fits 
have been made in the unpolarized case it is not possible yet in the polarized 
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case. Therefore we have convoluted the NLO polarized densities with the 
NNLO coefficient functions to get an idea of the stability of the perturbation 
series. First we have seen that the dominance of the sea-quark initiated pro- 
cess, which holds at LO and NLO, also prevails at NNLO. For the unpolarized 
cross section we have shown that the addition of the NNLO terms produces 
a significant reduction in the scale variation of the cross section. However 
the same phenomenon is also shown for the polarized cross section in spite 
of the fact that the NNLO parton densities are not known yet. Finally the 
longitudinal asymmetry is small and positive for the polarized parton density 
sets BB1, BB2 and GRSV(SS). If one uses the GRSV(VS) densities then the 
longitudinal polarization is negative and large. 
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Appendix A 



The lowest order contribution originating from the Born reaction Eq. (3.1) 
is given by 



*11 



-6(1 -z). 



(A.1) 



The 0(a s ) correction to the qq subprocess in Eq. (3.2) has been calculated 
in the literature [8]- [11]. Choosing the MS scheme the expression for A^ 
can be obtained by using n-dimensional regularization. For convenience we 
split the contribution as follows 



A&(z) 



AiF +V W+A^ H ( 



(A.2) 



Here S+V indicates that the corrections are coming from soft-plus-virtual 
gluons where H means that the corrections are coming from the hard z ^ 1 
gluon region. The expressions are equal to 



,(i),s+v 
Hi 



= -a,C F \S{\ 



6 In 



'Q 2 



+ 8C(2)-16 



+8V (z) \n(^ 



+ IQVAz) 



(A.3) 



,(1),H 
*qq 



= a s C F h{l + z) In (Qpj + 8(1 + z) ln(l - z) 



+4 \nz 

1 — z 



(A.4) 



where the distibutions T>i(z) are defined in Eq. (3.22). The second order 

(2) 

correction to the non-singlet part A q g is determined by process in Eq. (3.4) 
and the higher order corrections to processes in Eq. (3.1) and Eq. (3.2). It 
can be split into two parts. The first piece is related through mass factoriza- 
tion to the collinearly singular part of the partonic structure function W qq - 
and will be denoted by 



(2),NS _ a(2),S+V A (2),C A , a(2),C f »(2) a (2) 

11 r ^11 r 11 qq,AA ^ qq,AC 
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'■II 



A (2) - 

11, AD ■ 



(A.5) 
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(2) 

The second piece consists of the DY correction terms A _ d5 , A ( t'~,, A v _ ,,,, 

r qq,BB' qq,BC^ qq,BD 

(2) 

The contribution A - A § is zero because of Furry's theorem. First we will 
enumerate the contributions in Eq. (A. 5). The soft-plus- virtual gluon con- 
tributions come from the two-loop virtual graphs contributing to the Born 
reaction in Eq. (3.1), the soft gluon corrections to the reactions in Eqs. 
(3.2), (3.4) and soft quark pair production due to diagrams A in Fig. 1. The 
expression for this part is equal to 



(2) 
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+ 



192 V 2 (z) + 96 - (128 + 64 C(2)) V (z) 
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The hard gluon contribitions are denoted by A^' Ca and Afq'° F . They are 
equal to 
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+ yln 2 (z) - 121n(z) ln(l-z) 
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The function Li„(z) and S„ iP (z) denote the polylogarithms and can be found 
in [31]. The hard part of quark pair production due to the diagrams A in 
Fig. 1 is equal to 
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Finally, we have the interference terms corresponding to the combinations 
AC and AD in Fig. 1 and Fig. 2. For them we find 
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(A.10) 



The remaining parts of qq scattering are free of mass singularities. Therefore 
they do not need mass factorization, which implies that their contributions 
are scheme and scale independent. The contributions from the diagrams B 
in Fig. 1 and the interference terms BC and BD (see Figs. 1, 2) are 
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Because of Furry' s theorem 



i {2) - 



= 0. 



(A.13) 



At 0(a s ) the qg subprocess shows up for the first time. The Drell-Yan 
correction term for the reaction in Eq. (3.3) has been calculated in [8], [9], 
[10], [11] and it is given by 



Az-2 



ln ( (1 ^ 2g2 ) +(1 ~ z)(5 + 3z) }- (A ' 14) 



The second order part of A qg receives contributions from the virtual correc- 
tions to the reaction in Eq. (3.3) and the process in Eq. (3.5) 
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The calculation of requires both mass factorization and renormalization. 
The two parts A^ 2) ' Ca and A( 2),Cf are equal to 
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The reactions in Fig. 2 describe quark-anti-quark as well as quark-quark 
scattering (without identical quarks). The contribution to the DY correction 
term can be split into two parts. The first part, represented by the combi- 
nations CC and DD needs mass factorization. In this case the contributions 
for qq, qq and qq are all equal and are given by 






4(1 + z)\n(z) - 10(1 -z) In 2 




+ [8(1 + z) [ - 2Li 2 (l - z) + ln 2 (^) - 2 ln(z) ln(l - z) 



+4(7 - 8z) ln(z) - 40(1 - z) ln(l - z) + 40(1 - z) 



xln pU _40(l-^)[ln 2 (l-^)-C(2)] + 8(1 + 




x - 6Si )2 (l - z) + 4Li 3 (l -z)- ln(z)Li 2 (l - z) 




+2 ln 2 (^) ln(l - z) - 2 \n(z) ln 2 (l - z) - 4(9 - 7z) 



xLi 2 (l - z) + 8(7 - 8z) ln(z) ln(l - z) + (-25 + 39z) 



35 



x ln 2 (z) + 80(1 - z) ln(l - z) + 2(-42 + 17z) ln(z) 



-131 + 1283 + 3^ 



(A.18) 



The second part, which is collinearly finite, consists of the interference be- 
tween the graphs C and D in Fig. 2. Further notice the relative minus sign 
between the qq and the qq (qq) part. The expresions for these interference 
terms are 



A {2) - 

qqCD 



-A (2) - = -A (2) - 

qqCD qqCD 



a 2 C F T f 8(2 + z) - 4Si )2 (l - z) + 12S 1;2 (-z) 



+6 ln 2 (l + 2) ln(z) + 6C(2) ln(l + z) - 5 ln(l + z) ln 2 (z) 

+ 12 ln(l + z)U 2 (-z)] - 8(1 + z) [2Li 2 (-z) + 2 ln(*) ln(l + 2) 

+C(2)] - 8(2 - z)U z (l -z) + 16(6 - 5z)Li 3 (-z) 

16 

+8(6 - 9z)C(3) - 8(2 + 3z)Li 2 (l - ^) ln(^) + — ^ln 3 (^) 

o 

-128 ]n(z)U 2 (-z) - 8(2 + 5z)C(2) \n(z) - 8Li 2 (l - 2) 



+Az \n 2 (z) + 16 ln(z) + 32(1 - z) 



(A.19) 



In case there are identical quarks in the final state, we have in addition 
to the graphs in Fig. 2 also the ones in Fig. 3. As the results for EE, 
FF and EF are equal to those for CC , DD and CD (of course one has 
to implement the right statistical factors), we will not discuss them here. 
The new contributions come from the intereference terms CE, CF, DE and 
DF. Further when the intermediate vector boson is a photon all four sets 
of diagrams C,D,E and F contribute and we have a statistical factor 1/2. 
Apart from the statistical factor there is another difference bewtween CE 
{DF) and CF (DE). The first contains collinear divergences and needs mass 
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factorization, whereas the latter is free of mass singularities. The correction 
corresponding to the interferences CE and DF is equal to 



i (2) - 

qq,CE 



.(2) 



.(2) 



(2) 



qq,DF qq,CE qq,DF 



a 2 Cp [ Cf — t;Ca 



1 

-i 

2 



l + z 2 
l + z 



X 



4 ln 2 (^) - 8C(2) - 16Li 2 (-z) - 16 ln(z) ln(l + z) 



+8(1 + z) ln(z) + 16(1 - z) 



l + z 2 



ln {% ' l + z 



X 



32Si, 2 (l - z) - 16Si, 2 (-z) - 32Li 3 (l - z) 



-8Li 3 (-z) - 32Li 3 + 32Li 3 (i— ^) - 4((3) 

+24Li 2 (l - z) \n(z) + 32Li 2 (-z) \n(z) - 32Li 2 (-z) ln(l - z) 
-16Li 2 (-*) ln(l + z) + 12C(2) \n(z) - 16C(2) ln(l - z) 



-8C(2) ln(l + z) - - ln 3 (^) + 8 ln 2 (^) ln(l - z) 
o 

+28 \n 2 (z) ln(l + z) - 8 ln 2 (l + 2) ln(z) - 32 ln(z) ln(l - 2) 



x ln(l + z) 



+ (1-*) 



-16S 1)2 (-z) + 8Li 3 (-z) + 8C(3) 



-16Li 2 (-z) ln(l + z) + 4C(2) ln(z) - 8((2) ln(l + 2) 
4 ln 3 (^) + 4 ln 2 (^) ln(l + z) - 8 ln 2 (l + z) \n(z) 



+32m(l-z) -34 



+ (! + «) 



8Li 2 (-z) +4C(2) 
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+16 hv(z) ln(l - z) + 8 \n(z) ln(l + 2) 



+ 8(3 + 2) 



xLi 2 (l - z) - 4(1 + 3z) \n 2 (z) - 2(9 - 7z) ln(z) 



(A.20) 



The expression for the interference terms CF and DE is 



i {2) - 



.(2) 



,(2) 



qq,DE qq,CF 



1 (2) - 
qq,DE 



a 2 s C F {C F --C A )Ul-zf 



8Li,(l-z) 



-12Li 2 (l - 2) - 8Li 2 (l - 2) \n{z) - 8Si, 2 (l - 2) 
4 



--ln 3 (^) -61n 2 (^) 
3 



2(7-6*) ln(z) 



-(15 + 13^ 2 -28z) 



(A.21) 



Notice that the above expression is scheme and scale independent. The 
diagrams for the gluon-gluon subprocess can be obtained from the quark- 
anti-quark annihilation graphs via crossing. The subprocess shows up for 
the first time at 0(a 2 s ). We have divided its Drell-Yan correction term into 
two parts i.e. 



A(2) = a(2),c a + a(2),c f . 
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99 



<!<! 
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The C^-contribution is collinearly finite and is therefore scheme and scale 
independent. It is given by 



a(2),Ca 

99 



N 2 



N 2 - 1 



l + z) 2 16Si, 2 (-z) + 24Li 3 (-z) 



1 

+ 16 C(3) + — Li 2 (-z) - 24 Li 2 (-z) \n(z) 



+16Li 2 (-z) ln(l + z) + 8C(2) ln(l + z) + -((2) 

o 
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16 

-12 ln 2 {z) ln(l + z) + 8 ln 2 (l + z) ln(z) + — ln(z) 



x ln(l + z)] - 8(1 - z) 2 S h2 (l - Z ) + {A+ ^z + 22z 2 ) 



x 



\n(z) + ( 



4 2£ 



14 



3 3 Z 3 Z 



! )ln 2 (,) + ^(l-, 2 )}. (A.23) 



The CV-contribution contains collinear singularities. After mass factorization 
in the MS-scheme we find 



Ag>' c ' = a 2 s \\2(l + Az)\n(z) + 8(l-z) 



'-(¥, 



+ 



2(1 + 4s) (4 \n{z) ln(l - z) - \n 2 (z) + 4Li 2 (l - z)) 



+32(1 - z) ln(l - z) - 2(9 - Sz) ln(z) - 43 + 52z - 9z 2 
x In (^pj + 4(1 +Az)[- 4Li 3 (l - z) + Li 2 (l - z) \n{z) 

+4Li 2 (l - z) ln(l -z)- \n 2 (z) ln(l - z) + 2 ln(z) 
x ln 2 (l - *)] - 8(1 + z) [Li 2 (-z) + ln(z) ln(l + z) 
+4(1 + z) 2 [4Si, 2 (-z) + 4Li 2 (-z) ln(l + z) 
+2((2) ln(l + z) - 3 ln 2 (z) ln(l + 2) + 2 ln 2 (l + 2) 
x ln(z)] + 8(3 + 6z + ^ 2 )Si i2 (l - *) + 8(1 + 2z - z 2 ) 
xU 3 (-z) + 4(1 + 2z - 2z 2 )((3) + 4(3 - 4s)Li 2 (l - z) 
-8(2 + 4z + z 2 ) \n(z)U 2 (-z) - 4(3 + 10z + 2z 2 ) 



xC(2) Hz) - 4(9 - 7z)C(2) + 2(1 + -z + -z 2 ) ln 3 (z) 
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+2(7 + 3z) In 2 z + 32(1 - z) ln 2 (l - z) - 4(9 - Bz) 
x ln(z) ln(l - z) + (55 - 72,2 + 27^ 2 ) \n(z) 

-2(43 - 522 + 9z 2 ) ln(l - 2) + 2(35 - 31z - 4^ 2 ) I . (A.24) 
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Figure Captions 

Fig. 4. The NLO differential cross section dAa/dQ at y/s = 200 GeV for 
the BB1 parton densities plotted in the range 2 < Q < 30 GeV and 
fi 2 = Q 2 . The NLO plots are for the 0(a s ) corrected qq (dot-dashed 
line), the 0(a s ) contribution qg (dashed line) and the total sum (solid 
line). Note that all contributions have negative regions at large Q so 
we have plotted their absolute values (see text). 

Fig. 5. Same as Fig. 4 but now for NNLO. Further we have shown the 0(a 2 ) 
corrected qq (dot-dashed line), the 0(a 2 ) corrected qg (long-dashed 
line), the 0(a 2 ) contributions gg (dotted line) and qq (short-dashed 
line) and the total sum (solid line). Note that all contributions have 
negative regions so we have plotted their absolute values (see text). 

Fig. 6. The cross section dAa/dQ at ^/s = 200 GeV plotted in range 2 < 
Q < 30 GeV; LO (dotted line), NLO (dashed line), NNLO (solid line). 
Note that all contributions have negative regions at large Q so we have 
plotted their absolute values (see text). 

Fig. 7. The K-factor at = 200 GeV plotted in the region 2 < Q < 30 
GeV; K NLO (dashed line), K NNLO (solid line). 

Fig. 8a. The cross section dAa/dQ at \/s = 200 GeV plotted in the region 
2 < Q < 30 GeV; LO ii = Q (solid line), LO ix = Q/2 (dashed line) , 
LO n = 2Q (dotted line). 

Fig. 8b. The cross section dAa/dQ at yt~s = 200 GeV plotted in the region 
2 < Q < 30 GeV. NLO n = Q (solid line), NLO n = Q/2 (dashed 
line), NLO fi = 2Q (dotted line). 

Fig. 8c. The cross section dAa/dQ at \fs = 200 GeV plotted in the region 
2 < Q < 30 GeV. NNLO fi = Q (solid line), NNLO // = Q/2 (dashed 
line), NNLO fx = 2Q (dotted line). 

Fig. 9a. The unpolarized cross section da/dQ at yfs = 200 GeV plotted 
in the region 2 < Q < 30 GeV. LO ii = Q (solid line), LO ix = Q/2 
(dashed line), LO ix = 2Q (dotted line). 
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Fig. 9b. The unpolarized cross section da/dQ at y/s = 200 GeV plotted in 
the region 2 < Q < 30 GeV. NLO (i = Q (solid line), NLO fx = Q/2 
(dashed line), NLO fx = 2Q (dotted line). 

Fig. 9c. The unpolarized cross section da/dQ at y/s = 200 GeV plotted in 
the region 2 < Q < 30 GeV. NNLO fx = Q (solid line), NNLO fx = Q/2 
(dashed line), NNLO fx = 2Q (dotted line). 

Fig. 10a. The polarized quantity N(fx/fx ) (Eq. (4.4)) plotted in the range 
0.4 < fx/fx < 2 with Q — 5 GeV and fi = Q. The results are shown 
for LO (dotted line), NLO (dashed line), NNLO (solid line). 

Fig. 10b. Same as in Fig. 10a but now for Q — 10 GeV. 

Fig. 10c. Same as in Fig. 10a but now for Q — 15 GeV. 

Fig. 11a. The unpolarized quantity N(fx/fx ) (Eq. (4.4)) plotted in the 
range 0.4 < fx/fx < 2 with Q = 5 GeV and fi = Q. The results are 
shown for LO (dotted line), NLO (dashed line), NNLO (solid line). 

Fig. lib. Same as in Fig. 11a but now for Q = 10 GeV. 

Fig. 11c. Same as in Fig. 11a but now for Q — 15 GeV. 

Fig. 12. The longitudinal asymmetry A LL (Eq. (4.5)) in % plotted in the 
range 2 < Q < 30 GeV at fx = Q. LO (dotted line), NLO (dashed 
line), NNLO (solid line). 

Fig. 13. The longitudinal asymmetry All (Eq. (4.5)) in % plotted in the 
range 2 < Q < 30 GeV at fi = Q and NLO. BB1 (solid line), BB2 
(dot-dashed line), GRVS01, standard scenario (dashed line), GRSV01, 
valence scenario (dotted line). 

Fig. 14. The longitudinal asymmetry All (Eq. (4.5)) in % plotted in the 
range 2 < Q < 30 GeV at fx = Q and NNLO. BB1 (solid line), BB2 
(dot-dashed line), GRVS01, standard scenario (dashed line), GRSV01, 
valence scenario (dotted line). 
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